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Abstract 



If signals of new physics are discovered at the LHC it wih be crucial to determine the spin structure of the new model. We 
. . jdiscuss a method that can help to address this question with a low integrated luminosity, C = 1 fb~^, at ^/s = 14 TeV. 
Based on the differences in angular distributions of primarily produced particles we show that a significant difference can 
be observed in the final state jet-pairs rapidity distance. An additional advantage of the method is that it does not rely 
on any particular structure of the couplings in the decay chain. We simulate samples for models with supersymmetric 
^ and UED-like spin structure and show that a distinction can be made early on. 
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. The Large Hadron Collider (LHC) has started operat- 
^ r"! "ing and will soon probe physics at the TeV scale, perhaps 
Q^revealing the origins and mechanism of the electroweak 
Q^Tsymmetry breaking. One of the most promising candi- 
dates for explaining this phenomenon is supersymmetry 
r*| (SUSY) [l|, l2|. In supersymmetric theories each Standard 
Model (SM) particle is paired with a superpartner of spin 
different by 1/2. In particular, spin- 1/2 quarks will be ac- 
^ companied by spin-0 squarks, and spin-1 gluons by spin- 
'1/2 gluinos. Another possibility is provided by models 
.with universal extra dimensions (UED) [3]. In such mod- 
lels each SM particle will have a tower of different mass 
[Kaluza-Klein (KK) partners of the same spin. 

Since different models of new physics predict differ- 
ent spins for the newly discovered states, the determina- 
tion of spins will be of extreme importance for establishing 
„the new theory. Generically, SUSY and UED also have a 
different mass structure, but model-independent measure- 
inents of both masses and spins will be required in order 
to get handle on the underlying model. In this Letter 
we consider the possibility of distinguishing models with 
the same mass structure but the spin structure of either 
SUSY or UED. We will not refer here to any particular 
UED model, assuming only that it inherits all the proper- 
ties (masses and couplings) of the analyzed SUSY scenario 
apart from the spin. Whilst we are only interested in the 
generic spin structure, we will refer to the same-spin part- 
ners as KK-particles for ease of notation. As a benchmark 
we choose two mSUGRA derived scenarios. 

So far there have been a number of features suggested 
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that would hint at the particular spin structure at the LHC 
(see also [4]): 

• The total cross section [5^. 

• Observation of higher KK modes ^ . 

• Kinematic distributions of quarks from quark part- 
ners decays [tI, Isj . 

• Particle production in vector boson fusion 

• Invariant masses of lepton-jet [IS IHl E E Q 

and 

lepton-photon [iB^airs in squark/KK-quark and 
gluino/KK-gluon 16| decay chains. 



• Kinematic reconstruction of missing momentum pj, 

Angular distributions of leptons from sleptons [iqI , 
13, Hi I or s-channel resonances decays [22[, and b- 
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jets from bottom squarks decays [23]. 

Most of these methods require significant statistics and 
consequently a high luminosity (typically C ~ O(IOO) fb~^) 
or a very specific decay chain. Here, we propose an exten- 
sion of the method originally proposed in [19| for sleptons, 
to the first and second generation squarks/KK-quarks. In 
a large class of models, strongly-interacting states will pro- 
vide the first observation of any new physics. Therefore, 
they offer the opportunity to get hints of the spin structure 
with early data. Here, we consider an integrated luminos- 
ity of £ = 1 fb~^ at 14 TeV center-of-mass energy (cms), 
however, the method is also applicable at lower center-of- 
mass energies. 

The Letter is organized as follows. In the next section 
we define our observable and present a proof of concept 
that it might be useful in studying spin structure of the 
underlying model. In Section [3] we discuss our benchmark 
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Figure 1: Polar distributions of (a) squarks, (b) KK-quarks, Eq. (c) gluinos, (d) KK-gluons, Eq. (|2|), in the hard process CM frame 
normahzed to the respective total cross sections at = 14 TeV. Contributions from left and right states have been summed over. For 
simplicity we take rriq = rrig = 500 GeV, however, a departure from this assumption does not change the qualitative behaviour of the cross 
sections. Note that, as dictated by the hard process matrix element, the curves with distinguishable particles in the initial and final state are 
not symmetric with respect to the incoming parton but will be symmetric with respect to the proton beam. 



scenario and details of the event simulation. Section S] con- 
tains the results of the simulations and discussion. Finally, 
we conclude in Section [5l 

2. Spin-sensitive observable 

The production of gluinos and squarks of the 1st and 
2nd generation is in many scenarios a dominant source of 
super symmetric particles. It is, therefore, one of the most 
promising channels for SUSY searches at the LHC j24 |. 
We will show in this Letter that with early data (1 fb""*^) 
we can already deduce important hints about the spin of 
the produced particles. We focus here on the di-jet channel 
(i.e. at least two hard jets, see Sec. H] for details) for which 
the sample sources at the parton level are 



pp Qi q) 
pp^gg, 



pp^qiQj 



(1) 
(2) 



as well as the charge conjugated processes (for squarks), 
followed by the decays e.g. 



Qi QXn^ 



where i, j = L, i?, n = 1, . . . , 4, /c = 1, 2 and q' denotes a 
quark of different flavour. Of course, the decay chains will 
very much depend on the details of the spectrum of the 
model. One particular example leading to a di-jet signal 
common in mSUGRA scenarios is 



PP^ QrQr^ QQX^ Xi 



whereas in UED the respective process is 



PP qm qm ^ qq 71 71 



(4) 



(5) 



qqi 



etc.. 



(3) 



where qm and 71 are the KK-partners of the right-handed 
quark and the photon, respectively. However, many dif- 
ferent processes from Eqs. ([I]) + (|3]) and (j2j) + (|3j) can con- 
tribute to the di-jet final state and therefore, we consider 
a fully inclusive signal. In addition, at a hadron collider, 
extra QCD jets will appear when we include initial state 
radiation (ISR) and final state radiation (FSR). 

Depending on the spin of the primarily produced par- 
ticle, the distributions in the center-of-mass (CM) frame 
of the system will be substantially different. Here we fo- 
cus on the comparison between SUSY models and those 
with a UED-like spin structure where the quark partners 
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Figure 2: Normalized polar distributions of squarks (black solid) 
and KK-quarks (red dotted) in the hard process CM frame for rriq = 
rrig = 500 GeV with cms 14 TeV and with all contributions from 
strong-interacting states summed over. 



Figure 3: Parton level distribution of cos^*g, Eq. ((6)), for SUSY 
(black solid), Eq. (01), and UED (red dotted), Eq. in the pp CM 
frame for rriq = rrig = 500 GeV and m ~o = 100 GeV at = 14 TeV. 



are spin-0 and spin- 1/2, respectively, whereas the gluon 
partners carry spin- 1/2 and spin-1, respectively. Distri- 
butions for the production cross section, process ([1]) and 
(|2]), at the leading order have been calculated in the liter- 
ature and can be found in e.g. [25., 21, 28|, S |30| for 
SUSY and in for UED. Since the shape and relative 

size of the production channels depends on the collision 
energy, we fold in the relevant parton density functions 
(PDFs). This gives an effective polar CM angular distri- 
bution with respect to the beam at the LHC in both cases, 
see Fig.lH The distributions in Fig. [1] have been calculated 
at the leading order, but we note that for SUSY, the next- 
to- leading order (NLO) result is known [30]. The shape 
of the distributions remain nearly unchanged though and 
only the total cross sections differ. 

Different channels of squark/gluino (KK-quark/KK- 
gluon) production exhibit a different angular dependence 
(see e.g. [29]), however when comparing SUSY and UED 
one clearly sees that the production of supersymmetric 
particles (spin-0 quarks and spin- 1/2 gluinos) tends to 
be more central compared to KK-particles (spin- 1/2 KK- 
quarks, spin-1 KK-gluons), Fig.[Tl The difference between 
the two spin structures is shown clearly in Fig. [2] by sum- 
ming over all channels. The main contribution to this ef- 
fect comes from the fact that, depending on the final states 
under consideration, SUSY particles will be produced in 
either an S- or P-wave excitation. Particles produced in 
a P-wave, due to angular momentum conservation, will 
have the polar distribution vanishing for cosOq^ = ±1. 
On the other hand, KK-particles will be dominantly pro- 
duced in 5- waves where the distribution has a minimum at 
cosOqj^ =0, for more details see [32, 33]. This conclusion 
is independent of the particle masses as the KK-particles 
are produced in a more forward direction for all channels 
discussed. In addition, the produced particles will be heav- 
ily boosted and hence this will affect any decay products 
originating from the initial particle. 

To probe the production distribution more directly, we 



propose the following observable, originally suggested for 
slepton production [19|, but adapted here for squarks, 

cos (9*^ = tanh (^^y^^ ' ^^^^ " ~ . (6) 

where Ar]qq is the difference of the pseudorapidities be- 
tween the two final state quarks from squarks or KK- 
quarks decay chains, e.g. Eq. (jH) and (O, Fig. [3l As 
discussed in [19] this variable is the cosine of the polar 
angle of quarks with respect to the beam axis in the frame 
where the pseudorapidities of the quarks are equal and op- 
posite. Being a function of the difference of pseudorapidi- 
ties, it is longitudinally boost invariant. This observable 
has proven useful in studying the spin of sleptons [19] and 
sbottoms [23]. Here we apply a similar approach to the 
case of the 1st and 2nd generation squarks. 

The distribution of the quarks in the laboratory frame 
is also affected by the decay process. In the case of spar- 
ticle decay, the angular distribution of quarks will be flat 
as we have a decaying scalar. 



da 



oc 1 , 



(7) 



where cos^* is the quark polar angle with respect to the 
squark momentum, boosted to the squark rest frame. The 
situation is more complicated for KK-fermions. The distri- 
bution will now depend on both the chiral coupling struc- 
ture and the polarization of the KK-quark. For a polarized 
fermion we have 



da 



d cos 0* 



(X 1 ± 



fl 



fl 



(8) 



where Jl.r are left and right-chiral electroweak couplings, 
respectively. For = fn we retain the distribution from 
Eq. ([7]) and there is no distinction from a scalar decay. For 
the other extreme case, /l = or /^^ = 0, the distribution 
will have a triangular shape. However, since fermions (KK- 
quarks in this case) will not be fully polarized, in general 
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the net distribution of quarks will be ex 1 ± a cos 6* with 
< a < 1. Another issue is that anti-fermions will have 
an oppositely shaped distribution, cf. ± sign in Eq. ([8]). 

Even if the KK-quarks have a high degree of polariza- 
tion, it should be noted that at the LHC, the particles are 
produced with a significant boost. Therefore the dominant 
structure of the distributions will be due to the produc- 
tion angle of the initially produced particles. Any angular 
structure from the decaying particles can only be expected 
to change the distribution by a small amount regardless of 
the coupling structure. In contrast, the majority of meth- 
ods so far proposed for spin determination at the LHC rely 
on a specific electroweak coupling structure in the model 
(namely Jl ^ fn or Jl ^ /i?) [l^- They would also re- 
quire correct identification of particles in the specific decay 
chains. 



3. Simulation and benchmark points 

In order to reahstically assess the observability at the 
LHC, we perform a hadronic level simulation using Her- 
wig++ 2.4.2 [si, [H, [36| with MRST 2004L0 PDFs [13 • 
Events are then analyzed using Rivet 38, [sol 
anti-A:^ jet algorithm ^ |4l| with Ai? = 0.5. 



and the 
The jets 

are required to be in the central part of the detector, i.e. 
\r]\ < 2.5. 

We include the following SM backgrounds in the anal- 
ysis: VK^+jets, Z+jets and ti. Matrix elements are gener- 
ated using MadGraph [42] and hadronized using Herwig++. 
We do not include the and di-boson backgrounds 

here as these are found to be subdominant [24] for the di- 
jet signal after the cuts listed in Sec. HI As benchmark 
scenarios we choose the ATLAS SU3 and SU6 mSUGRA 
parameter points 0. The SU3 scenario features lower 
masses and is slightly above the current ATLAS and CMS 
exclusion limits [43|,|4J, with squark masses of 0(650) GeV. 
For the slightly heavier SU6 scenario the squark masses 
are in the range of 0(850) GeV. Masses and parameters 
for both scenarios are given in Tab. [H The total NLO 
SUSY cross sections are 22.9 pb and 6.2 pb for SU3 and 
SU6, respectively. We simulate the samples of signal and 
backgrounds corresponding to an integrated luminosity of 
1 fb""*^ at ^/s = 14 TeV for both scenarios. Due to differ- 
ent spins, the total cross section for UED is larger than for 
SUSY assuming the same masses (by a factor ^ 10 in our 
scenarios). Therefore, we normalize the number of UED 
events to the number of SUSY events at NLO calculated 
using Prospino [30] for our comparison. 



^QCD background is subject to a large theoretical uncertainty 
0(100%). The correct treatment would require full detector simu- 
lation, inclusion of higher-order effects and finally the real data. If 
it turns out that the QCD background is too large, one can use the 
2j + i channel, but this in general would require a higher integrated 
luminosity [2^ . 
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Table 1: mSUGRA parameters and particle masses (in GeV) 
for SU3 and SU6 scenarios 24]. Masses were calculated using 
SPheno 2.3 [45j. 



4. Numerical results 

In our analysis we focus on the 2-jet signal for which 
we require at least two hard jets but no hard leptons in 
the final state. In order to suppress the QCD background, 
relatively high-p^ jets are required along with a hard EJ^^^^ 
cut. For each of the scenarios we employ the following set 
of cuts 

• At least two jets with p-l^ > 150 GeV and > 
100 GeV for the hardest and second-hardest jet, re- 
spectively. 

• No electrons and muons with pt > 20 GeV and 
\Vi\ <2.5. 

• A(j){ji,E^''') > 0.2, A^(j2,^T''') > 0-2. 

• E^'^^ > O.SMeff. 

• Meff > 800 GeV, 

with the effective mass Meff defined as 

M,^=p^ ^p^' ^Ef^' . 

A detailed ATLAS study Q , including a detector sim- 
ulation, showed that the above mentioned set of cuts give a 
good signal-to-background ratio (~ 2 — 8) for SUSY mod- 
els with ruq < 900 GeV and high statistics even at low 
integrated luminosity ^ 1 fb~^. Therefore, in the present 
study we focus on the possibility of seeing hints of the spin 
structure of the new physics at 14 TeV LHC in the early 
data. However, the same method can be applied at lower 
center-of-mass energies, albeit with lower statistics. 

Our observable, Eq. ([6]), has to be correspondingly 
modified at the jet level to be 



cos Ojj = tanh 



(9) 



where ji and j2 are the hardest and 2nd hardest jets, re- 
spectively. We are analyzing an inclusive SUSY/UED sig- 
nal, therefore all the channels from Eqs. ([T]) and ([2]) will 
contribute. However, this turns out not to be a problem as 
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Figure 4: Distribution of pseudorapidity difference for jets, Eq. (|9l), in SU3 ((a) and (c)) and SU6 ((b) and (d)) scenarios for SUSY (black) 
and UED (red) spin structure. The dotted blue line is the SM contribution. The number of events corresponds to an integrated luminosity 
C = l fb-^ ((a) and (b)) and £ = 10 fb-^ ((c) and (d)) at 14 TeV cms for SUSY particle production. The number of events for UED has 
been normalized to the number of SUSY events, respectively. The l-cr statistical error on each point is shown. 



the gluon partners production and the associated quark- 
gluon partners production ah possess the desired features 
observed for process @ and (j5j), see Fig. [TJc) and[Hd). 

Note that after the inclusion of hadronization and QCD 
radiation we are no longer sure whether the observed hard 
jet originates from a parton from the hard process. This 
is not a problem in the SUSY-like scenarios discussed here 
because of the large mass hierarchy, but may be a diffi- 
culty in other types of models. For example, if we take 
the simplest UED model, the mass splitting between the 
KK-quarks and the KK-photon is normally much smaller 
than in the models presented here [46]. Consequently, the 
jets produced when the KK-quarks decay will typically be 
relatively soft. Isolating these jets will therefore be signif- 
icantly more difficult and will certainly require a different 
set of cuts than those used for this study. 

Figure H] shows that even for simulated events, large 
differences remain between a model with SUSY and UED- 
like spin structures for both scenarios, left and right panel, 
respectively. In Fig. ID^a) and Hl^b) we expose the clear 
differences present with 1 fb~^ of data. Figures ID^c) and 
md) show how the situation becomes even clearer after 
10 fb"^ have been collected. In Fig. [5] we show the limiting 
case of very high statistics for SUSY and UED, where the 



qualitative difference between the models is clearly visible. 

Several differences can be noted when comparing the 
hadron-level. Figs. H] and O and the parton-level distri- 
butions. Fig. [3l However, these differences can be easily 
understood once the effect of the experimental cuts have 
been considered. The first difference is that the number of 
events where cosOjj ~ ±1 is lower than the parton-level 
expectation for both SUSY and UED. This is the result of 
a rapidity cut |7^| < 2.5 on the final state jets. Thus, large 
rapidity differences between jets are less likely as one or 
both jets in the observable will not be reconstructed. 

We also notice a dip in the distribution centered around 
cosOjj ~ 0. The dip is due to the fact that the recon- 
structed jets have a finite size. Two jets with a small 
rapidity difference can therefore overlap if they are close 
in the azimuthal direction and will not be resolved as be- 
ing separate. Thus, events with this topology will not be 
reconstructed. 

To quantify in one number a clearly visible difference 
between the distributions we introduce the following asym- 
metry, 

7V(|cos^*.| >0.5)-A/-(|cos^*.| <0.5) ^ ^ 
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Figure 5: Comparison of distributions between SUSY (solid black) 
and UED (dotted red) for the limiting case of very high luminosity 
for the point SU3 at ^/s = U TeV. 



where N{. . .) is the number of events fulfihing the re- 
spective condition. For our benchmark scenario SU3 at 



1 fb we obtain, 



/lobs 



^SUSY 

/lobs 
■*^UED 



-0.22 ±0.01 



= 0.01 ± 0.01 



SUSY — -0.226 , 
UED= 0.016, 



where ^susy' ^ued high statistics hmits of the 

asymmetry for SUSY and UED, respectively, and ^susy' 
-^UED observed values with statistical errors. Thus, 

using the outlined method we can clearly distinguish be- 
tween the two spin structures investigated. 

Finally, we consider different possible mass hierarchies 
of squarks and gluinos (KK-quarks and KK-gluons). The 
original SU3 and SU6 scenarios feature similar masses of 
quark and gluon partners, see Tab.[Tl Here we also analyze 
two extreme cases, with either rriq <C rrig or rriq rrig. 
In Tab. [2] we have collected the values of the asymmetry, 
Eq. (p!Q|) . for different mass ratios of quark partners to 
gluon partners. All other parameters are the same as in 
the SU6 scenario. The values of the asymmetry ^ ^ 
and A ~ —0.4 would clearly point to a particular scenario. 
However, if ^ ^ —0.2 we are left with an ambiguity, as this 
value can be obtained in both UED and SUSY scenarios. 
In principle, the final spin determination would require at 
least an approximate knowledge of the mass scales in the 
new physics sector. 

Nevertheless, the case with heavy squarks/KK-quarks, 
rriq ^ m^, is special since the di-jet signal is produced 
here by 3-body decays of gluinos/KK-gluons. In particu- 
lar, neither the cuts nor the observable itself are optimized 
for this type of the signal. Clearly, in order to come to a 
decisive conclusion about the spin, an additional observ- 
able would have to be included to confirm or exclude 3- 
body decay case. For example, if gauginos decay to leptons 
a clear invariant mass, m^^, edge would point to 2-body 
quark partner decays. On the other hand, 3-body decays 
of gluon partners would contain 4 partons in the final state 
with the similar pt distributions. In contrast, with quark 
partners decays, only two high-p^ partons occur in the fi- 
nal state. Therefore, triggering on additional jet activity. 





rriq <C rrig 


rriq ^ rrig 


rriq > rrig 


UED 


0.05 





-0.19 


SUSY 


-0.22 


-0.22 


-0.38 



Table 2: Comparison of asymmetries, Eq. (|10|) , in SUSY and UED 
for different mass scenarios of quark and gluon partners. Other pa- 
rameters are the same as in the SU6 scenario. The statistical error is 
0.02, assuming the number of events as in the original SU6 scenario 
and an integrated luminosity C = 1 fb~^ . 



together with the study of angular distributions and invari- 
ant masses, should in principle allow us to separate these 
two cases. We also note that heavy quarks/gluon partners 
would significantly lower the cross sections, hence more 
integrated luminosity would be required. 

5. Conclusions 

We have presented a method that may help to ex- 
pose the spin structure of new physics that may be seen 
in early LHC data. By comparing the distributions of 
the di-jet rapidity distance for the models of SUSY and 
UED spin structure, we found that this observable may 
provide good supporting evidence to discriminate between 
them. The method relies only on the spin of initially pro- 
duced strongly-interacting particles. It does not require 
any particular electroweak coupling structure and it works 
for many possible decay chains. Therefore, it can be ap- 
plied to a wide range of scenarios. The results of the sim- 
ulation show a clear distinction between models already 
with C = l fb~^ Although we have performed the Monte 
Carlo simulation for ^/s = 14 TeV the same conclusions 
hold at lower center-of-mass energies, albeit with lower 
statistics. While more elaborate studies will be needed in 
order to fully confirm the spins of all individual particles 
in a model, early hints of the spin structure will help to 
direct further measurements. 
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